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TOWARDS OPTIMAL REGULARITY FOR THE 
FOURTH-ORDER THIN FILM EQUATION IN 
GRAVELEAU-TYPE FOCUSING SELF-SIMILARITY 

P. ALVAREZ-CAUDEVILLA, J.D. EVANS, AND V.A. GALAKTIONOV 


Abstract. An approach to some “optimal” (more precisely, non-improvable) regularity of so¬ 
lutions of the thin film equation 

ut = —V • (|w|"VAu) in x R+, M(a;, 0) = rio(®) in R”^, 

where n £ (0,2) is a fixed exponent, with smooth compactly supported initial data uo{x), 
in dimensions Y > 2 is discussed. Namely, a precise exponent for the Holder continuity with 
respect to the spatial radial variable |x| is obtained by construction of a Graveleau-type focusing 
self-similar solution. As a consequence, optimal regularity of the gradient Vu in certain 
spaces, as well as a Holder continuity property of solutions with respect to x and t, are derived, 
which cannot be obtained by classic standard methods of integral identities-inequalities. Several 
profiles for the solutions in the cases n = 0 and n > 0 are also plotted. 

In general, we claim that, even for arbitrarily small n > 0 and positive analytic initial data 
uo(x), the solutions u{x,t) cannot be better than C'|“^-smooth, where e(n) = 0{n) as n —>■ 0. 


1. Introduction: TFE~4 and known related regularity results 

This paper is devoted to an “optimal regularity” analysis for a fourth-order quasilinear degener¬ 
ate evolution equation of the parabolic type, called the thin film equation (TFE-4), with a fixed 
exponent n > 0, 

(1.1) ut = —V • (juPVAu) in x M+, u{x,t)) = uq{x) in , 

where n G (0, 2) is a fixed exponent, and bounded, sufficiently smooth, and compactly supported 
initial data uq with an arbitrary dimension N >2. 

Moreover, we assume uq = no(|x|) to be radially symmetric, so solutions u = u{\x\,t) do the 
same. These initial conditions could be relaxed (for example uq G (1 , etc.). However, for 

simplicity, we have chosen those initial conditions, since the focusing phenomenon to be studied 
here exists in any reasonable functional class of initial data, including radial ones. 

The main result obtained here is an approach to the so-called optimal regularity of solutions of 
the TFE-4 for the A-dimensional case with N >2, ascertaining an “optimal” Holder continuity 
exponent in x and t. 
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For the one-dimensional case, similar results were obtained in [7], establishing the Holder 
continuity in C ^'2 with respect to the variable x and in with respect to t. However, in 
the A^-dimensional case, the optimal regularity has been unsolved besides the interest of the 
specialised mathematical community; see li nni la US], in many of those works it is also 
assumed that the solutions are non-negative for the Cauchy problem (CP). However, as we 


recently proved in the CP settings, the solutions for the TFE~4 (1.1) are oscillatory and sing¬ 


changing for all not that large values of n > 0; see 121 El for recent results and references 
therein. 


One of the difficulties concerning optimal regularity for the equation (1.1) is due to the fact 
that this equation (1.1) is not fully divergent, so the whole “global” regularity theory for the 


TFE-4 can use just a couple of well-known integral identities/inequalities. Indeed, as discussed 
in [2| an integral identities argument, as that performed by Bernis-Eriedman in ID, fails to 
show the Holder continuity in higher dimensions; see [2| for details and reasons for these issues. 
Consequently, we claim that, for such “partially” (and/or fully non-) divergent equations, a 
different approach should be put in charge. 

1.1. Methodology: Graveleau-type focusing similarity solutions. Here, working on the 
radially symmetric problem of TFE-4, such that 

u{x,t) = u^{r,t), with r = \x\, 

we base our analysis on the foeusing argument performed, firstly, by Graveleau [Hj in a formal 
sense, and later justified in Aronson-Graveleau [I] for a class of non-negative solutions (via the 
Maximum Principle, MP) of the elassic porous medium equation (PME-2): 

(1.2) vt = A(u™') in x (—oo,0), with an exponent m > 1. 


Actually, in comparison with (1.1), m = 1 -|- n > 1. In particular, the authors in [3] obtained, 
after constructing a one-parameter family of self-similar solutions, that, 

when N >2, there exists 0 < ^ = g.{m,N) < 1, 

and a radial self-similar solution u* (r, t) to the focusing problem and with the focusing time 
t = 0~, such tha10 

(1.3) u*(r,0“) = C^r^, where C* is an arbitrary constant. 

This hnished a long discussion concerning optimal regularity for the PME-2 after seminar results 
by Caffarelli-Friedman [9| (see MathSciNet for many further extensions) proving, earlier. Holder 
continuity of v{x,t) in a very general setting. 

The above results of Aronson and Graveleau showed that such a Holder continuity is, indeed, 


optimal, and even the corresponding Holder exponent p, = p{m, N) from (1.3) can be evaluated. 


at least, numerically for any m > 1 and N > 2. Indeed, the blow-up singularity (1.3) shows. 


actually, not an optimal regularity for such PME-2 flows, but the non-improvable one, in the 


sense that, for the Holder exponent of solutions of (1.2), 

(1.4) a non-improvable regularity =► Holder exponent p. 

Fortunately, as we have mentioned above, Caffarelli-Friedman earlier results of the 1980s had 


been already available at that time, so, together with (1.4), those proved that a proper Holder 


regularity of any non-negative solutions of the PME-2 is optimal. 


^We do not put further details how to get (1.31, since we will explain this shortly for the TFE-4. 
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These Aronson-Graveleau focusing (singularity-blow-up) ideas for the PME~2 (1.2) and other 
second-order quasilinear parabolic equations (e.g. for the p-Laplacian one) got further devel¬ 
opment and extensions in a number of papers devoted to various important aspects of porous 
medium flows; see references and results in 13 El- 

More precisely talking about such Graveleau-type similarity approach, from the point of view 
of applications, in the focusing problem, one assumes that, initially, there is a non-empty compact 
set K in the complement of the suppuo, where v vanishes. In other words, there is a hole in 
the support of the initial value u(x,0) = t'o(x) > 0, and, in finite time T, this hole disappears, 
making the solution v to become positive along the boundary of K and eventually at all points 
of K. Basically, due to the finite propagation property that these equations possess (the porous 
medium equation and the thin film equation). Thus, as the flow evolves the liquid enters K and 
eventually reaches all points of K at the instant T. We are then interested in how the solution 


for the TFE-4 (1.1) behaves near the focusing time T. We again suppose (as in the Graveleau’s 
argument for PME-2) the focusing time t = 0~. 

However, in our case, with the TEE-4 (1.1) and n > 0, but close to zero, we need to take into 


consideration the existence of oscillatory solutions of changing sign. 


Furthermore, for the TFE-4 (1.1), there are not still any general regularity results even in the 
radial setting in . Consequently, we will follow the lines and the logic of (1.4), i.e. we will 


estimate a certain non-improvable Holder exponent for radial (and, hence, all other) solutions. 


To this end, we use those previous singularity (blow-up) ideas from the PME-2 (1.2) to 
establish some properties for the self-similar solutions of the radially symmetric problem of 
the TFE-4 (1.1) and, eventually, some optimal regularity information about its more general 


solutions. Indeed, as we will see through the analysis performed in this paper, to obtain such 
an “optimal” (non-improvable) regularity for the TFE-4 (1.1), it seems that we must solve a 
nonlinear focusing problem for that TFE-4 that will be derived from the associated self-similar 
equation. 

To be precise, we first apply the focusing ideas performed by Aronson-Graveleau [1] for the 
PME-2 (1.2) to the TFE-4 (1.1). To do so, we will work on the space of radially self-similar 
solutions with r = |x| > 0 of the form 


(1.5) 


y = 


r 

(±IF 


for ± t > 0, where /3 = 


l-\-an 


> 0 . 


Thus, these solutions will solve an associated self-similar equation, or nonlinear eigenvalue prob¬ 
lem, 

-J=T[y^-Vr{J=T{y^-^ry)y± ^yf ^af = o, 

with — or -|- depending on if we are analysing the focusing phenomena or the defocusing phe¬ 
nomena. In other words, before the focusing time or after (in our situation t = 0). 

As usual, the unknown a priori a represents the nonlinear eigenvalues supporting the reason 
to call that equation a nonlinear eigenvalue problem. 

Crucial in our analysis will be to impose some radiation-type conditions, or some minimal 
and maximal growth at infinity, i.e. as y —)• -|-oo. This actually allows us, among other things, 
to determine the existence of a discrete family of nonlinear eigenvalues denoted by {afc}- 

Moreover, this behaviour at infinity leads us to have the solutions of the self-similar TFE-4 
bounded at infinity by functions of the form 


f{y) = Cy^{l+o{l)) C G M and /x =/x(a, n) := > 0. 
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Thus, we ascertain a limit at the focusing point 

u^{r,t) ^ Cr>^ as t— >-0“, with /x = ^, 

uniformly on compact subsets (see details in Section]^. However, this does not provide yet 
information about the influence of the dimension on the behaviour of the solutions at infinity. 

To get to that point, we need to find those nonlinear eigenvalue^ Wk} associated with a 
discrete family of nonlinear eigenfunctions {fk}- Here, we have done this analysis via a homotopy 
deformation as re —)• 0^ (Sections and to the linear problem with re = 0, 

(1.6) ut = —A^re in x M_, 


and with those patterns occurring at re = 0 acting as branching points of nonlinear eigenfunctions 
of the Cauchy problem (1.1) when re is close to zero. In particular, according to our analysis, 
when re = 0, it follows that 

aA:(0) = Q!fc = |, k = 1 , 2 , 3 , - . 

As an observation, those values of the parameter a when re = 0 might be written as a perturba¬ 
tion of the eigenvalues for the fourth-order operator 

B*f = -\yf = A/, 

analysed in full detail in m-, see further comments below. 

Hence, we perform a homotopy deformation from the TFE-4 (1.1) to the parabolic bi¬ 
harmonic equation (1.6), for which we can apply again a similar logic in ascertaining Graveleau- 
type “focusing solutions”. Indeed, we get a minimal and maximal growth for the radial self¬ 
similar solutions of this linear problem (1.6). Moreover, since we know that there exists a discrete 
family of eigenvalues for the corresponding self-similar equation associated with (1.6), using this 
branching/homotopy argument we obtain the desired family of values for the parameter a from 
which we will have a family of radial self-similar solutions emanating at re = 0 from the self¬ 
similar solutions of the linear problem (1.6). Note also that the family of eigenfunctions for the 
linear problem (1.6) is a complete set of generalised Hermite polynomials with finite oscillatory 
properties HU. 

In this paper, we also support this analysis numerically, performing a shooting procedure from 
y = 0 to y = -|-oo that gives us those linear eigenvalues ak and the profile of the corresponding 
eigenfunctions. 

Additionally, we show some numerical analysis that provided us, with very difficult to ascer¬ 
tain, profiles of the nonlinear eigenfunctions {fk} with re > 0 and sufficiently close to zero. 

This numerical analysis supports the conjecture that by homotopy continuity the properties 
for the profile when re = 0 remain valid for re > 0 and sufficiently close to zero. As seen above, 
we have fixed the interval re G (0,2), however, the numerics suggest that those properties are 
lost from re = 0.8 or 0.9. In fact, if re > 2 we just have nonnegative solutions as shown in [7] 
so that, this focusing argument is not possible. Thus, we restrict this analysis to re G (0,2), or 
more precisely, of re > 0 sufficiently close to zero. 


^Since we obtain those nonlinear eigenvalues as a perturbation from a family of linear eigenvalues when the 
parameter n is close to zero, we will use indistinctively this notation for the family of linear or nonlinear eigenvalues. 
To be precise, the nonlinear eigenvalues correspond to the family {ak{n)} and the linear eigenvalues to the family 
{Qfc(O)}. 
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Finally, we obtain such an “optimal” (non-improvable) regularity for the TFE-4 (1.1) proving 
that there is Holder continuity with an specific coefficient 


^^k — 


4afc 


1 + nak 


and Vu(x, t) G with p<p*{n,N) = 


N 


1 


having such regularisation depending on the dimension N, as we claimed and were looking for. 
Consequently, if n > 0 small, we find that 


(1.7) 


t{r,t) £ with £{n) = npL + 0{n^). 


It turns out that, thanks to our branching analysis, there holds /x > 0, so that, for n > 0, 
the regularity condition becomes obviously worse. Note that, if /U < 0, we have a better 
regularity than for n = 0 which is impossible. 

In general, we obtain (and the numerics suggest that as well) that there exist several focusing 
singularities in the radial geometry of the type with A: = 1, 2,3, • • •. In particular, for the 

values of ak obtained here, it follows that for Ac = 1, ai = ^ is the minimal and crucial one, 
since it seems to be the only degenerate case, i.e. /i(0) = 0. All others satisfy 

/fc(0) >0, A' = 1, 2, 3, • • • , i.e. the TFE-4 for n > 0 is not degenerate 


initially, but, indeed, will be eventually, at the focusing time t = T~{= 0“). Therefore, we can 
conclude that slight changes in the parameter a destroy the regularity, supporting the fact that 
these results are valid only when n is sufficiently close to zero. Indeed, via branching analysis 
we find that 

ak{n) = ak + Pkn + O(n^), 
so that we have the regularity condition 

U,{r,t) £ C‘^k-i,,n+Oin^) ^ 


Note that, even for a 2 = 1, the analytic positive solution becomes, at t = T , i.e. not 

classical solutions in Hence, for the minimal ai = ^ it is which is much worse. 


Remark. As mentioned above, we also observe that, for the PME-2 (1.2), one needs to have 
such a non-empty compact hole to apply the Maximum Principle. However, we believe that, for 
the TEE~4 (1.1) and the value ai = g, that provides us with the minimal regularity 
there is not such a hole. We do not have a rigorous justification of it but the numerics presented 
in this paper suggest it that way. 


In relation to the Holder regularity with respect to the temporal variable t, assuming the 
radially self-similar solutions of the form (1.5) we find that 

|n(0,t) -m( 0,0)| = (-t)“/(0). 


Hence,the Holder’s exponent for the variable t when it is very close to t = 0 cannot be bigger 
than 

ai{n) = 1 / 2 - 1 - 0{n). 

Therefore, we provide with an estimation from above for the Holder continuity with respect to 
the t variable. Moreover, this estimation improves the Holder’s exponent obtained by Bernis- 
Eriedman [7] showing that theirs was actually not optimal since it was 1/8. 
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1.2. TFE-4 problem settings. This setting is well-known nowadays (though many things were 
not fully proved in a general case), so, below, we omit many details; see surveys [T5l|T6]. Note 
that, for some values of n > 0 (not that large), focusing similarity solutions to be constructed 
do not exhibit finite interfaces, so the resulting optimal regularity results are true for any FBP 
and/or Cauchy problem settings. Principal differences between the CP and the standard FBP 
settings for the TFE-4 (1.1) are explained in nails]. 

We recall that the solutions are assumed to satisfy the following zero contact angle boundary 
conditions: 


( 1 . 8 ) 


u = 0, 

Vu = 0, 

—n • V(|tt|” Au) = 0, 


zero-height, 

zero contact angle, 

conservation of mass (zero-flux) 


at the singularity surface (interface) TqIu], which is the lateral boundary of 


supp u C 




N > 1 


where n stands for the unit outward normal to Tq [u] , which is assumed to be sufficiently smooth 
(the treatment of such hypotheses is not any goal of this paper). For smooth interfaces, the 
condition on the flux can be read as 


Next, we denote by 
(1.9) 


lim —n • V(|n|"Au) = 0. 

dist(a:,ro[n])10 


M{t) := f u(x, t) dx 


the mass of the solution, where integration is performed over the support. Then, differentiating 
M[t) with respect to t and applying the divergence theorem, we have that 

= - I n-V(lurAu). 

ron{t} 


The mass is conserved if J(t) = 0, which is satisfied by the flux condition in (1.8) 


2. Radial self-similar solutions: focusing and defocusing cases 
Now, we construct the operators and specific solutions in order to apply the ideas performed by 


Aronson-Graveleau [1] for the PME-2 (1.2) to the TFE~4 (1.1). 

Thus, thanks to the scaling invariant property of these nonlinear parabolic equations, we now 
construct radially self-similar solutions, i.e. in terms of r = |x| > 0, with a still unknown value 
of the parameter a > 0 (clearly, it must be positive, as shown below) 


( 2 . 1 ) 


'-(r,t) = (±t)‘^f(y), 2 /=^^ for ± t > 0, where 13 = 


- ^ > 0 . 


Here, by the time-translation, we ascribe the blow-up or focusing time to T = 0 . We then 
simultaneously consider two cases: 

(i) Focusing (i.e., Graveleau-type) similarity solutions, which play a key role, corresponding 
to (—t) in (2.1) and the singular blow-up limit as t —>■ 0“, and 


11 


Defocusing similarity solutions, with {+t) in (2.1), playing a secondary role as exten¬ 


sions of the previous ones for t > 0, i.e., corresponding to the not-that-singular (or, at 
least, less singular) limit t —)• O'*'. Actually, these defocusing solutions are well-posed 










solutions of the Cauchy problem for the TFE-4 for t > 0 with initial data u* (r, 0 ) 
obtained from the previous blow-up limit as t —)• 0 “. 


Substituting (2.1) into (1.1) we arrive at the similarity profiles f{y) satisfying the following 
nonlinear eigenvalue problems: 


(2.2) B±(a, /) = -Vy ■ {\f\-VyAyf) ±(3yVyf^af = 0 for y > 0; /'(O) = /'"(0) = 0, 

where Vy and Ay stand for the radial gradient and the radial Laplacian. In ( |2.2| ), we present 
two symmetry conditions at the origin (which should be modified if f{y) = 0 near the origin, 
i.e. it contains a “zero hole” nearby). As we will see with the numerical analysis performed in 
Section of this paper, we can choose other conditions depending on if either / = 0 or / 7 ^ 0. 

To complete these nonlinear eigenvalue settings one needs extra “radiation-type” (or growth- 
type) conditions at infinity to be introduced next. 

Indeed, we actually find radially similarity profiles / depending on the single variable y. The 


operator of the equation ( 2 . 2 ) is then 


(2.3) 


A±(a, /) = -J^[y^-Vr{J^iy^-^fy)y± Pyf T «/ = 0, 


denoting the radial nonlinear operator as 

(2.4) R„(a,/) = ^[y^-vr{^{y^-^fyyt 

Thus, here, a > 0 is a parameter, which, in fact, stands in both cases for admitted real 
nonlinear eigenvalues to be determined, in the focusing Graveleau-type case, via the solvability 
of the corresponding nonlinear eigenvalue problem, accomplished with some special (radiation- 
type) conditions at infinity, to be introduced shortly. 


Remark 2.1. In general, with respect to the similarity profiles described above we note that 
there are two main types of self-similar solutions. For solutions of the first kind the similarity 
variable y can be determined a priori from dimensional considerations and conservation laws, 
such as the conservation of mass (1.9) or momentum. 

For solutions of the second kind the exponent /3 (and by relations the exponent a) in the 
similarity variable must be obtained along with the solution by solving a nonlinear eigenvalue 
problem of the form ( 2 . 2 ). 

The first published examples of self-similar solutions of second kind are due to G. Guder- 
ley in 1942 m studying imploding shock waves, although the term was introduced by Ya. 
B. Zel’dovich in 1956 [l 8 j. Further examples might be found in the theory of the collapse of 
bubbles in compressible fluids or in works on gas motion under an impulsive load; see Barenblatt 
[ 6 ] for an extensive work on this matter. 


3. Minimal growth at infinity (a “nonlinear radiation condition”) for 

Graveleau-type profiles 


Here, we consider the blow-up problem (2.3)_, i.e. with the lower signs. One concludes that, 
in order to get a possible discrete set of eigenvalues a > 0 , some extra conditions denoted by 
nonlinear radiation condition on the behaviour (in particular, growth) of f{y) as ?/ —)• -|-oo must 
be imposed. 
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Obviously, such a “radiation-type” condition (we use a standard term from dispersion theory) 
just follows from analysing all possible types of such a behaviour, which can be admitted by the 


ODE (2.3)_, which is not that a difficult problem. There are two cases: 

(I) More difficult: there is a “zero hole” for f(y) near the origin. Finite interfaces for the 
TFE-4 are well known in the case of the standard FBP setting and also in the Cauchy 
one (see dans] and references therein). Then we need to look for profiles /(y) which 
vanish at finite y = yo > 0 and describe asymptotics of the general solution satisfying 
zero contact angle and zero flux conditions as y ^ y^ 


f{y)^0, f{y)^0, -|/r/"'(y) 


0 . 


This case causes a difficult problem, since the oscillatory behaviour of f{y) close to the 
interface is quite tricky for the CP da, while, for the FBP, it is much better understood. 
(II) More standard and easier: /(O) 0. This case also includes the border possibility 

yo = 0, i.e. when /(O) = 0, but / ^ 0 in any arbitrarily small neighbourhood of y = 0. 

3.1. Minimal and maximal growth at infinity. This is key for our regularity analysis. Our 
radial ODE (2.3)_ admits two kinds of asymptotic behaviour at infinity, i.e. when yo —)• -|-cx). 


For the operator A„ {a, /), we state the following result: 


Proposition 3.1. For any a> 0, the ODE (2.3)_, with the operator A„(a,/) possesses: 
(i) Solutions f{y) with a minimal growth 

(3.1) /(y) = Cy^(l-ko(l)) as y ^+oo, CgR, 

where 


Aa 


> 0 . 


(3.2) y = y(a, n) = 

(ii) Moreover, there exist solutions of ( |2.3[ )_ with a maximal growth 

(3.3) /(y)~y''°, as y ^+ 00 , where /^o = |- 

Remark. 


It is important to mention that the expression (3.3) for solutions of the maximal growth 
does not include the “additional or extra” corresponding oscillatory component ip{s), 
with s = Iny, and shows just the growth behaviour of its “envelope” with algebraic 
growth. Such an oscillatory maximal behaviour will be introduced and studied later. 
Note also that, as n —?• O'*', we have that 


(3.4) 


4 

n 


+C)0, 


which corresponds to an exponential oscillatory growth for the linear problem occurring 
for n = 0; see next section. 


The difference between (3.1) and (3.3), which implies such terms as minimal and maximal 


(3.5) 


growth at infinity, is obvious: 

for any a > 0, 




4a 

l+an ‘ 


Proof. This follows from a balancing of linear and nonlinear operators in this ODE, though 
a rigorous justihcation is rather involved and technical. A formal derivation is surely standard 
and easy: 
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(i) In this first case, we assume a linear asymptotic as y —)• +oo (assuming simple radial 
behaviour / 

(3.6) i^y/'-a/ + ... = 0 ^ /(?/) ~ where ^ = ^>0. 

Since this behaviour is asymptotically linear, eventually, we get an arbitrary constant 


C 7^ 0 in (3.1). A full justification of existence of such orbits is straightforward, since the 
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nonlinear term in the ODE is then negligible. So that in the equivalent integral equation 
though a singular term, it produces a negligible perturbation. 

On the other hand, the solutions are also bounded by a maximal growth, which in this 
case comes from non-linear asymptotics that balance all three operators 


(3.7) 


R.n(/) + ^2//'-«/ = 0 


f{y) ~ y" as y ^ +00, 


where we again indicate the envelope behaviour of this oscillatory bundle. A justification 
here, via Banach’s contraction principle is even easier, but rather technical, so we omit 
details. □ 

Overall this allows us to formulate such a condition at infinity, which now takes a clear 
’‘minimal nature” such that solutions / are now bounded at inhnity by a function 


(3.8) 


f{y) = Cy^^{l + o[l)), with 


M = 4a ^ Q 
y l+an ^ 


Obviously, thus we just need a global solution f[y) of our ODE (2.3)_ in ]R_|_, satisfying the 


minimal growth (3.1). Indeed, for such prohles, there exists a finite limit at the focusing (blow¬ 
up) point: 

(3.9) u*(r, t) —>■ as t —>■ 0“, with 


t —>■ u , witn y = ^ 

uniformly on compact intervals in r = |x| > 0. One can see that, for any maximal profile as in 


(3.3), the limit as in (3.9) is infinite, so that such similarity solutions do not leave a finite trace 
as t —)■ 0~. 

However, the above proposition leaves aside the principle question on the dimensions of the 
corresponding minimal and maximal bundles as y —)• -|-cx), which is not a straightforward prob¬ 
lem. Note that the latter is actually supposed to determine the strategy of a well posed shooting 
of possible solutions of the above focusing problem. 

Indeed, the main problem is how to find those admitted values of nonlinear eigenvalues 
{oik}k>i (possibly and hopefully, a discrete set), for which / = fk{y) exist producing finite 
limits as in (3.9). 


To clarify those issues, we consider the much simpler linear case when n = 0 and, subsequently, 
pass to the limit as n —?• 0+ in (2.3). This analysis will provide us, eventually, with some 


qualitative information about the solutions, at least when n is very close to zero. 


4. The linear problem: discretisation of a for n = 0 


(4.1) 


For n = 0, the TFE-4 (1.1) becomes the classic bi-harmonic equation 

N 


ut = —Am 


m 


X M_. 


Of course, solutions of (4.1) are analytic in both x and t, so any focusing for it makes no 


sense. However, we will show that (following a similar philosophy to the one above) Graveleau- 


self-similar solutions of (1.1), at least, for small n > 0. 
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type “focusing solutions” for (4.1) are rather helpful to predict some properties of true blow-up 























4.1. Maximal and minimal bundles. Thus, first we consider the same “focusing” solutions 

(2.1) _ for (4.1), that take a simpler form 

(4.2) u*(r,f) = (-!)"/(?/), y = = 

Then, the corresponding linear radial ODE ( |2.3[ ) takes also a simpler form 

(4.3) /) = - s y/' + «/ = 0- 

with a maximal behaviour. Those are exponentially 


We first calculate the solutions of 
growing solutions of the form 


(4.4) 


f{y) ~ e“*^^ as y —>■ +oo 




This characteristic equation gives two roots with Re(.) > 0: 

(4.5) ai ,2 = I 4“i [4 ± i ^] = Co ± ici. 

and one negative root (this actually goes to the bundle of minimal solutions; see below) 

(4.6) 


as = -f 4 3 < 0. 


+00 (including the multiplicative 


Hence, the bundle of maximal solutions is oscillatory as y 
algebraic factor), 

(4.7) f{y) ~ y“i(^+^“)e'=o^'‘^®[Cicos(ciyt) + (72 sin(ciyt)], Ci,2 € M. 

In fact, obviously, this bundle is 3D since it includes the ID sub-bundle of exponentially decaying 
solutions with the exponent (4.6). However, in a shooting procedure performed at the end of 
this section, we intend to get rid of just two coefficients in (4.7) 

Ci = C2 = 0 . 

Indeed, this numerical evidence shows that those coefficients vanish at certain values of the 
parameter a that we will determine below specifically. Moreover, we provide an explanation in 
the next subsection. 


(4.8) 


On the other hand, the minimal behaviour (3.1) now reads 

/(y) = + 0(1 )), (7gM (with |=4a). 


The whole bundle of such minimal solutions is 2D. Besides the parameter C in (4.8), it includes 
a ID sub-bundle of exponentially decaying solutions with the exponent (4.6), so that, overall, 
the minimal solutions compose a 2D family: 

(4.9) /(y)~(7y4"(l + o(l)) + De“3?^"^'(l + o(l)), C, D G M. 

Justification of both behaviours is, indeed, simpler for such a standard linear ODE problem. 

Einally, passing to the limit as in (3.9) then gives, for minimal solutions, a finite “focusing 
trace”: 


(4.10) 


u* (r, t) —)• as t —> 0 


We explain why those trivial results are so important in what follows. The conclusion from 

(4.10) is straightforward: since n*(r, 0~) must be analytic then we have the following. 
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Proposition 4.1. For the above linear “focusing eigenvalue problem” (4.3), there exists not 
more than a countable set of admissible eigenvalues {afc}fc>i given by 


(4.11) 


Ofc = |, k = 1,2,3,.... 


Proof. The function is analytic at r = 0 only for values ak from (4.11), i.e. 4a must be 
a real even integer, say 2k. Note also that, in the present problem, all functions are analytic, 
so that, automatically, the set of roots {ofc} is discrete with a possible accumulation point at 
infinity only. □ 

Remark. Solving the linear eigenvalue problem: 


Ao,y(afc,/fc) = 0 in 


fk&Ll 


then, it seems that the nonlinear eigenvalue problem (2.3)_ formally reduces to the classic linear 
eigenvalue problem (4.3) at n = 0, providing us with another reason to call (2.3) a nonlinear 
eigenvalue problem. 

Also, the values of the parameter a can be written as shifting from the eigenvalues of the 
eigenvalue problem 

Wf = -\yf = A/, 

analysed in HH and whose discrete spectrum takes the form 


a{W) = {\k = -\-,k = D,!,■■■}. 


Hence, 


ak = \k-^ 


4 ’ 


with A; = 1, 2, 3, 


having a countable family of eigenvalues for the problem (4.3). 

Moreover, note that Aq y{a,f) (4.3) is a non-symmetric linear operator, which is bounded 


from Hp 


to Tp(M) with the exponential weight given by (4.7), i.e., 

p{y) = ^ Co > 0 small and defined by (4.5). 


4.2. Well-posed shooting procedure. We now discuss a practical procedure to obtain the 
linear eigenvalues a^. We perform standard shooting from y = 0 to y = +00 for the ODE (4.3) 
by posing four conditions at the origin: 


• Either 


(4.12) /(O) = 1 (normalisation), f''{ll) = 0, /^(O) = f'''{ll) = 0 (symmetry); 

• or 

(4.13) /(O) = 0, f''{ll) = 1 (normalisation), f'{ll) = = 0 (symmetry). 

These two sets of conditions correspond to partitioning of the eigenfunctions into two subsets 
with the corresponding properties at the origin of 

either /(O) = 1, /"(O) = 0, or /(O) = 0, /"(O) = 1, 

11 














together with symmetry conditions. Explicitly, the first four eigenfunctions are 


k = 

1 : 

CX\ — 

1 

2 ’ 

hiv) 

k = 

2 : 

0:2 = 

1, 

f2{y) 

k = 

3 : 

0:3 = 

3 

2 ’ 

fsiy) 

k = 

4 : 

a4 = 

2, 

h{y) 


where the above properties are immediately apparent. 

Numerically we may illustrate the appearance of the eigenfunctions and their eigenvalues 


through consideration of an Initial Value Problem (IVP). The linear ODE (4.3) is solved nu¬ 
merically using odel5s (with tight error tolerances AbsTol=RelTol=10“^^) and subject to either 


(4.12) or (4.13) as initial conditions. The far-held behaviour (4.7) is extracted from the numeri¬ 


cal solution. The oscillatory component in the far-held behaviour is revealed by considering the 
scaled solution 

(4.14) 


A least squares htting of this function to the remaining oscillatory component in (4.7) over a 


suitable interval for large y, allows the determination of the constants Ci^ 2 - The large y interval 
taken was typically [250,300]. The scaled function (4.14) is shown in Figure[^for each of the two 
types of initial conditions in the V = 1, 2 and 3 cases. Two selected values of the parameter a 
are taken in each case. The hgures illustrate convergence to the oscillatory part of the far-held 
behaviour, with the extracted least squares estimates of the constants Ci ^2 shown inset for each 
of the two a values. Figure shows the variation of the far-held constants Ci ^2 with a in the 
V = 1, 2, 3 cases and the two initial conditions. The constants C \^2 both vanish precisely at the 
eigenvalues, when a = Ok, the hrst hve being shown in Figure]^ in each N case. The magnitude 
of the constants Ci ,2 grow rapidly as o. increases, making determination of further eigenvalues 
more difficult. 

This approach demonstrates the validity of numerical determination of the eigenvalues and 


eigenfunctions by choosing a to minimise the far-held behaviour (4.7) 


5. A “homotopic” transition to small n > 0: SOME KEY ISSUES 

Next, we are going to use the above linear results to predict true nonlinear eigenvalues ak{n) 
for the TFE-4, at least, for sufficiently small n > 0. By continuity (to be discussed later on as 
a continuous deformation after applying a homotopic argument) we now know that 

(5.1) afc(0) = |, /j = 1,2,3,.... 

Another important conclusion; since, for n = 0, any /^(O) ^ 0 (this could happen only acciden¬ 
tally, with a probability 0), we can also expect that 

(5.2) for small n > 0, /a:(0) / 0, 

meaning that, in this case, we do not need to perform a shooting from the interface point, with 
a quite tricky behaviour nearby. 


Also, we conclude from (4.7) that: 

(5.3) for small n > 0, maximal bundle is 2D and oscillatory as y 

Those properties will aid progress on the nonlinear eigenvalue problem. 
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Figure 1. Numerical illustration of the oscillatory component in the far-held be¬ 
haviour (4.7) in the linear case. In each dimensional case N = 1,2,3 scaled prohles 


(4.14) are shown for two selected values of the parameter a for each initial condition 
(4.12) and (4.13). The extracted least squares values of the far-held constants Ci ^2 are 


stated inset in each hgure. 


5.1. Passing to the limit n —y O'*" in the nonlinear eigenvalue problem. Subsequently, we 


perform a homotopy deformation from the self-similar equation (2.3)_ (with lower signs) of the 


TFE-4 (1.1) to the linear radial ODE (4.3) corresponding to the classical bi-harmonic parabolic 


equation (4.1). In particular we construct a continuous deformation from the radial equation 


(2.3)_ to the linear equation (4.3) for which we know solutions explicitly. 
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(A) N=1,f(0)=0 (B) N=1,f(0)=1 






(E) N=3, f(0)=0 



(F) N=3, f(0)=1 



Figure 2. Numerical determination of the far-held constants Ci ^2 for (4.7) in the 
stated dimensional and initial condition cases. Shown are their variation with the pa¬ 
rameter a. The coincident zeros correspond to the vanishing of the far-held maximal 
bundle yielding the eigenvalues of the linear problem a = ak- 


It is clear that the CP for the bi-harmonic equation (4.1) is well-posed and has a unique 
solution given by the convolution 


(5.4) 


u{x,t) = b{x — ■,t) * uo{-), 


where b{x,t) is the fundamental solution of the operator Dt + A^. 

Due to our analysis is possible to establish a connection between the radial solutions of 0 
and the corresponding to (4.1) via the self-similar associated equations when n —)• O"*". To this end 
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we apply the Lyapunov-Schmidt method to ascertain qualitative properties of the self-similar 


equation (2.3)_ following a similar analysis as the one carried out in [T]. 

Thus, as we already know, the operator A, 
family of eigenvalues 

ak = afc(O) = |, with /c = 0,1, 




{a, f) defined by (4.3) produces a countable 


Note also that, (4.3) admits a complete and closed set of eigenfunctions being generalised Hermite 


polynomials, which exhibit finite oscillatory properties. 

This oscillatory issue seems to be crucial. In fact, in [T^] was observe that a similar analysis 
of blow-up patterns for a TFE-4 like (1.1) did not detect any stable oscillatory behaviour of 
solutions near the interfaces of the radially symmetric associated equation. Hence, all the 
blow-up patterns turned out to be nonnegative, which is a specific feature of the PDE under 
consideration therein. However, this does not mean that blow-up similarity solutions of the 
CP do not change sign near the interfaces or inside the support. Actually, it was pointed out 
that local sign-preserving property could be attributed only to the blow-up ODE and not to the 


whole PDE (1.1). Hence, the possibility of having oscillatory solutions cannot be ruled out for 
every case. 

5.2. Branching/bifurcation analysis. Now, we construct a continuous deformation such that 


the patterns occurring for the nonlinear eigenvalue problem (2.3)_ are homotopically connected 


to the ones of the equation (4.3). Thus, we assume for small n > 0 in (2.3) the following 
expansions: 

(5.5) ak{n) := ak +fJ-i,kn +o{n), |/|"-= l/l" = := 1-h nln |/|-b o(n), 

where the last one is assumed to be understood in a weak sense. The second expansion cannot 
be interpreted pointwise for oscillatory changing sign solutions f{y), though now these functions 
are assumed to have finite number of zero surfaces (as the generalised Hermite polynomials for 


n = 0 do). Indeed, as discussed in [T] for (2.3) this is true if the zeros are transversal. 


Furthermore, in order to apply the Lyapunov-Schmidt branching analysis we suppose the 
expansion 

(5.6) / = E|/3|=fc C/3//3 + Vk, for every A: > 1, 

under the natural “normalising” constraint 


E cg = l. 

m=k 


(5.7) 

Moreover, we write 

{//3}|/9|=fc {/l) •••) /Mj,}) 

as the natural basis of the M^-dimensional eigenspace, with Mk > 1, such that 

fk = E|/ 3 |=fc cgfg and 14 = E|/3|>fe C/3//3, 

with 14 G Tfc where Y). is the complementary invariant subspace of corresponding kernel. In 
particular, we consider the expansion 

14 := -b o(n). 

Thus, applying the Fredholm alternative, we obtain the existence of a number of branches 
emanating from the solutions {oik-, fk) at the value of the parameter n = 0. We can guarantee 
that the first profile is unique but for the rest there could be more that branch of solutions 
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emanating at n = 0. The number of branches will depend on the dimension of the eigenspace so 
that the dimension is somehow involved; see [T] for any further comments and a detailed analysis 
of a similar branching analysis. 

Remark. Furthermore, when n —)• 0"*“ we have a few different profiles of f{y). For the PME-2 


The previous discussion can be summarised in the following Lemma. 


( 1 . 2 ) 


( 1 . 1 ) 


the Graveleau profiles are always unique by the Maximum Principle but not for the TFE-4 
. However, for a = h the only profile is f{y) = y‘^. 


Lemma 5.1. The patterns oceurring for equation (4.3), i.e. the ones for (2.3)_ for n = 0, aet 
as branching points of nonlinear eigenfunetions of the Cauchy problem for the operator (2.3)_, 
at least when the parameter n is suffieiently close to zero. 


Remark. It turns out that using classical branching theory “nonlinear eigenfunctions” /(y) of 
changing sign, which satisfies the nonlinear eigenvalue problem ( |2.3| )_ (with an extra “radiation- 
minimal-like” condition at infinity) at least, for sufficiently small n > 0, can be connected with 
eigenfunctions fk of the linear problem (4.3). 


6. Non-improvable regularity for the TFE-4 
The following main result is a straightforward consequence of our focusing self-similar analysis. 


Theorem 6.1. Let, for a fixed n > 0, the nonlinear eigenvalue problem (2.3)_ have a nontrivial 
solution (eigenfunction) fk{y) for some eigenvalue > 0, i.e. there exists a self-similar focusing 
solution of the problem (2.2)_ exhibiting the finite-time trace (3.9). Then: 


(i) For the general Cauchy problem for the TFE-4 (1.1)? even in the radial setting, the 
Holder continuity exponent of solutions eannot exceec^ 


( 6 . 1 ) 




Tk{n,N) - pi - i+nofe- 


(ii) Let there exist < 1? be. au < for n G (0, 2)Fl Then, for the TFE-4 (1.1); 


( 6 . 2 ) 

(6.3) 


Vrt*(r,0 ) G Lf, 


loc 


iff pe[l,P*), P*{n,N) = 
so that, for (even radial) solutions of (O), in general, for any t > 0, 


Vn(x,t)0Lf;jM^). 


Remark. Basically |Vu*(r, 0 )| G Lf^^{M.^) if and only if 

/kv |Vn*(r, 0-)|P = Cfg r^-i(rf^k-ijP < qo. 

Thus, since 0 < y* < 1 we find that to have that integral bounded we need that, evaluating the 
coefficients in the second integral 

p<p,=p,{n,N) = 13 ^ 4 ^. 

^We mean if fj-k > 1, that, not that surprisingly, happens for all small n > 0. 

^Actually, for smaller n’s; Note that for every larger n’s solutions of the TFE-4 are known to be strictly 
positive, Cl). 
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Remark. Of course, it can happen that /r*, > 1, so focusing does not supply us with a truly 
Holder continuous focusing trace but rather than a C^+^-trace. Actually, by continuity in n, 
exactly this happens for n > 0 small, where fii{0,N) = 2. 

We expect that the minimal value of in (6.1) is attained at /c = 1, but cannot prove this for 
all n > 0 (for small ones, this is obvious). Note again that, for all large n > 2, such a focusing 
is not possible in principle, since the solutions must remain strictly positive for all times, [7]. 

Thus, this optimal (non-improvable) regularity results for the TFE-4 (and, seems for many 
other parabolic equations) depends on the solvability of the nonlinear eigenfunction focusing 
problem. 

Note that in the previous section we showed a homotopy connection between that nonlinear 
eigenvalue problem and the linear problem at n = 0 (4.3). For those reasons the analysis 
performed in those two previous sections is crucial in ascertaining qualitative information about 
the solutions of the TFE-4. 

Furthermore, we can also ascertain the Holder continuity with respect to the temporal variable 
t. This fact comes directly from the non-improvable Holder’s exponent we have already obtained 
above. Indeed, assuming the radial self-similar solutions of the form (1.5) for the focusing, one 
easily finds that 

|u(0, t) - m( 0, 0)1 = (-t)"/(0). 


Hence, the t-Holder’s exponent close to t = 0 cannot be bigger than a. 

Therefore, the focusing solution given here provides us with optimal Holder estimations for 
the variables x and t. 


7. Oscillatory structure of maximal solutions for n > 0 


Let us now consider n G (0,2). For the ODE (2.3), we will try the same anzatz as in |13j . 
Namely, we will find solutions of the maximal type, with the envelope as in (3.3). We introduce 
a corresponding oscillatory component as follows: 


(7.1) 

Then, 


fiy)=y^^is), s = lny, 


fa = — for y 1. 
n 


f = yt^ ^{(p + fj,ip), f" = y>^ {fa-l)ip + fi{fa-l)ip], 

f'" = ^ Yp + 3(^ — 1)(^ + (3^2 _ 0^ _|_ 2'^(p -\- _ 3 ^ _)_ 2) , 

where ' = d/dy and ‘ = d/ds. Substituting these expressions into (2.3) yields the following 
fourth-order homogeneous ODE for p>{s) 


(7.2) 

-h 2{N -4 + Ifi)'^ -h {6fi‘^ + 6{N - 4)fi + 11 + {N - 1){N - 9))(p 
-|- 2(2// -|- — 4)(//^ -|- (iV — 4)/i -|- 2 — N'jip -|- /i(/t — 2)(/i^ -|- 2[N — 3)// -|- 3 -|- [N — 1){N — 5'))ip 

— 4 -|- 3/i)(/3 -|- (3/i^ -|- 2[N — 4)/t -|- 4 — 2N)ip fa{fa — 2')[N — 2 -|- fa)<p\ 

-7 (1(1 -h na){ip + flip) - ap) \p\-'^ = 0. 


We mention that this equation is autonomous and thus may be reduced to third-order, although 
we do not utilise this reduction here. Figure [^illustrates the periodic nature of the solutions for 
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(f), at least for n small enough. Since the oscillations occur over such a large range, we use the 
following transformation to allow the oscillations to be visible on the plots, 


(7.3) 


t((?i(s)) = 


ln(?!)(s) + 1, if 4>{s) > 1, 

4>{s), if - 1 < (j){s) < 1, 

— ln(—(/>(s)) — 1, if (?i(s) < —1. 


Plotted are numerical solutions for cj) in the two cases a = 0.5,1 with N = 1. The subplots 
illustrate the change in the profile behaviour with n. The numerics support the proposition that 
the family (7.1) is 2D, composed of: 

(i) a ID stable manifold 
and 

(ii) a phase shift + sq) for any sq G K. 

As we have seen earlier, these two properties are true for re = 0, and, hence, by continuity, we 
conjecture remain true for small re > 0. 

However, the periodic exponential structure, with the linear ODE for re = 0 is replaced by a 
more difficult one (7.1) for re > 0. Nevertheless, the periodic nature of such a behaviour appears 
universal, and is expected to remain for larger re. The numerics though do suggest that it is lost 
(in a homoclinic-heteroclinic bifurcation) for re near to 0.8 or 0.9. 

To reconcile with the known behaviour in the linear case, we now study the behaviour of the 
periodic solutions for small re > 0. To reveal the limiting oscillatory behaviour as re —t- 0, we 
keep the leading terms of the coefficients in (7.2) as —>• oo to obtain 

(7.4) 


1 


^ + 4// 7? -I- if + f if + -{f + ff) \ f\ "■ = 0 


f{s) = f ^fis), 


s 

s= 


We now rescale as follows 

(7.5) 
leading to 

(7.6) f +4:lf+%lf + 4:lf + f+\{^if + i^\f\~'^ = Q, 


where here ’ denotes d/ds. For re = 0, the equation (7.6) becomes linear, 

(7.7) if + Af + Qif + Af + if + \ {f + =0, 

with the characteristic equation 

{m 1) {{m + 1)^ + i) = 0, 

for exponential solutions if{s) = e™'*. The roots of the characteristic equation are linked to those 
of the controlling factor in (4.4) being 


4 

3 ® 2 , 




m -|- 1 = |ai. 

Consequently we obtain the dominant asymptotic behaviour 


and 0. 


(7.8) 


f{s) 


,(-l+|co)s 


Ai cos (|ci.s) -|- A 2 cos (|cis) 


as s 


00 , 
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N=1,a=0.5 


(A) n=0.01 (B) n=0.1 




(C) n=0.5 


(D) n=0.8 





Figure 3. Illustrative numerical profiles of transformed 0(s) via (7.3). Shown is the 
parameter case = l,a = 0.5 for selected n. 


for arbitrary constants Ai^ 2 - Denoting the independent variable by y rather than y for conve¬ 
nience, we thus have the small n behaviour 

f{y) ~ v" (t)"" ^ i^Hy)) ■ 

19 


(7.9) 



















































N=1,a=1 


(A) n=0.01 


(B) n=0.1 



(C) n=0.5 



(E) n=0.9 




(F) n=1.5 



Figure 4. Illustrative numerical profiles of transformed (j){s) via (7.3). Shown is the 
parameter case = l,a = 1 for selected n. 


This may be reconciled with the expression in (4.7) through the identifications 

M„4/3 -- ... 3 




8. Nonlinear eigenfunctions by shooting 


The eigenfunctions for n > 0 may be obtained via shooting. The conditions (4.12) or (4.13) are 


again used as initial conditions and a determined by capturing the growth (3.8) for sufficiently 
large y values. This growth behaviour may be imposed by 

minimising fiyf' — af at y values, 
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typically chosen to be around 40. Standard regularisation of the term |/(y)|” is required in the 
form (/^ + (5^) 2 with 6 taken relatively small. 

Figure shows the eigenvalues of the hrst three branches k = 1,2,3. These numerically 
extend the n = 0 eigenvalues of Proposition |4.1| to n > 0. 

As in the n = 0 linear case, the eigenvalues remain the same irrespective of the spatial 
dimension N, at least to the accuracy of the numerical calculations that were performed. 

Furthermore, it is worth remarking that along each branch i.e. fixed k, the exponent fi of the 
far-held behaviour of the eigenfunctions 

/ ~ Cy^, 

remains quite close to its value in the linear case i.e., 

fj.(n) « /i(0). 

Thus, as an approximation we have 

ak(n) « afc(O)/(1 - afc(O)n), 
which seems to be a reasonable approximation at least 

for n < l/ak(0). 

The corresponding eigenfunctions are thus similar to those in the linear case. 


N=1 N=2 




Figure 5. Numerical determination of the eigenvalues akin) for the first three 
branches k = 1, 2,3 of eigenfunctions. 


An alternative approach to determining the eigenfunctions in both the linear n = 0 and 
nonlinear n > 0 cases, is through minimisation of the oscillatory maximal profile to obtain the 
non-oscillatory minimal profile for y ^ 1. 

In principle we have two parameters a and, say v = f"{Q) or /(O), depending on which branch 
of eigenfunctions is considered. In the linear case, using (4.7), we are required to satisfy two 
algebraic equations with analytic functions: 


(Ciia,v) = 0 , 

\C'2(a, u) = 0. 
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Therefore, we arrive at a well-posed 2 — 2 shooting problem, which cannot have more than a 
countable pairs of solutions (as mentioned already). This approach is practical for the linear 
n = 0 case, as the two parameter form of the maximal bundle is known explicitly as given in 
(4.7) and was essentially pursued in section 4.2. However, lacking such an explicit expression 
for n > 0, limits this approach in the nonlinear case. 

Finally it is worth presenting some numerical experiments showing the behaviour of the 
maximal profiles. Figure shows profiles in the case N = l,a = 0.75 for selected n. Since the 
oscillatory profiles occur over such large ranges, the transformed profile using (7.3) is depicted. 

The figures suggest the loss of a pure oscillatory structure as n increases, which is compatible 
with the behaviour seen for (p in Figures and 

We conjecture that this is highly suggestive of a homoclinic-heteroclinic bifurcation. The 
precise values of n though at which this occurs is difficult to determine with sufficient accuracy. 


9. After-focusing self-similar extension 


As usual and as we have observed in many similar problems on a (unique) extension of a 
solution after blow-up, this is much easier. In fact, after focusing at f = 0“, we arrive at a 
well-posed CP (or TFE, that does not matter since uo{r) > 0, no interfaces are available) for 
the TFE-4 0 with initial data ( |3.9[ ) (with fi = when possible), already satisfying the 
necessary minimal growth at infinity. 

Therefore, there exists a self-similar solution (2.1)+ with f{y) satisfying the ODE (2.3)+, with 
already fixed “eigenvalue” a = a^. Because of changing the signs in front of the linear terms, 
this changes the dimension of the stable/unstable manifolds as y —>■ -|-oo (in particular, the 
unstable manifold becomes ID, similar to n = 0). 

To explain the latter, as usual consider the case n = 0 (and hence small re > 0). Then, with 
the change of sign in the two linear terms in (4.3), we arrive at a different characteristic equation 
for a: 


(9.1) /( 2 /)~e“S'^ as y ^+oo 

This characteristic equation gives two stable roots with Re(-) < 0: 


(9.2) fti,2 = f 4 3 [ — 4 ± i — Co ± ici. 

and one positive root maximal solutions; see below) 

63 = I 4"5 > 0. 


Hence, the bundle of maximal solutions is oscillatory as y —)• -|-oo is just ID, so we arrive at 
an undetermined problem for re, a, which satisfy just a single algebraic equation (unlike two in 
. Indeed, this makes the a-spectrum continuous and solvability for any a > 0. 

This allows us to get a unique extension profile /(y) for such a a priori fixed eigenvalue ak¬ 
in other words, the focusing extension problem is not an eigenvalue one, since a proper has 
been fixed by the previous focusing blow-up evolution. 

It might be said, using standard terminology from linear operator theory, that for the sign -|- 
the spectrum of this problem becomes continuous, unlike the discrete one for —. Therefore, we 
do not study this much easier problem anymore, especially, since it has nothing to do with our 
goal: to detect an non-improvable regularity for the TEE“4 via blow-up focusing. 
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N=1,a=0.75 


(A) n=0 (B) n=0.01 




y y 


(C) n=0.5 (D) n=0.8 




{E)n=0.9 (F)n=1.5 




Figure 6. Illustrative numerical profiles of transformed (/)(s) via (7.3). Shown is the 
parameter case = l,a = 1 for selected n. 
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Appendix A. The limit n —f 0 for maximal solutions 


as 


We consider here the singular limit re —)> O"*" for the equation in (3.7), written explicitly here 


+ 2(1 + an)yf - af = 0. 


The non-uniform solution in this limit comprises two regions, an Inner region where 1 <C y < 
0(re“^/'^) with I/I” ~ 1, and an Outer region y = 0{n where In |/| = 0(1/re). The labelling 
of these regions as inner and outer becomes apparent during the course of the scalings. 

In the inner region 1 <C y < 0(re“^/^), we obtain at leading order in re the linear ODE (4.3). 
Posing / ~ /o with |/o|” ~ 1, we obtain 


(A.l) 




-h 1(1 + are)y/^ - afo = 0. 


The far-field behaviour of (A.l) may be determined using a WKBJ expansion in the form 
(A.2) 

which gives 


/o(y) ~ A(y)e‘^(^) as y -boo 


(A.3) 




1 

= 4 ^’ 


3yA' + 2 (A + 2a - 1 - 12(/)')^</>") ^ = 0. 
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The required solutions to (A.3) take the form 


(A.4) 


(p{y) = ay 






A{y) = ky 3 


where a satisfies the cubic equation in (4.4) and k an arbitrary constant. Thus, the dominant 
behaviour for large y is 

(A.5) fo{y) ~ (^/siexpjaiyt} +/c 2 exp { 022 / 5 }^ , 

with 01^2 as given in ( |4.5[ ) and ki ^2 complex constants chosen so that the expression is real as 
stated in (4.7). 

This inner solution breaks down when y = where ln|/o| = 0(l/n). This suggests 

3 

the consideration of an outer region with scaling Y = n~iy. In T = 0(1), we have 


(A.6) 


yiV-l dV 


yW-l|^|n d 


1 


dY V dV 


vAf-1 d/ 
^ dV 


+ \{l + an)Y^ - af = 0. 


Rather than posing a multiple-scales ansatz directly, it is more convenient to work in complex 
form and we may instead consider 


(A.7) 


fiY) 


.biy)/n 


B{Y) 


as n 


0 , 


where h is complex in order to match with the inner solution expression (A.5). Thus, at 0(l/n) 
in (A.6) we obtain 

(A.8) |e^| {b'f + = 0, 

whilst at 0(1) we have 

(A.9) 3f+ (l-f+ ln|R|)6' + 6f^ + ^i^M±^=0, 

where ' denotes ^ and the approximation 


(A.IO) 


i/r 


(1 + n In |i?|) 


has been used. The solution to ( A.8[ ) that matches with ( |A.2[ ) and (A.4) is 
(A.ll) 5(y) = oA In (1 + . 

The amplitude B{Y) is determined from (A.9) using the solution (|A.11[). We obtain 


InlRl = 


(i-r^W/3) 


- |(A + 2a)lny- g(2iV + 3a-4)y4/3) 


after matching to (A.2) with (A.4). 
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